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ABSTRACT. It is known that the Barlow-Proschan index of a system with i.i.d. component lifetimes coincides with the 
Shapley value, a concept introduced earlier in cooperative game theory. Due to a result by Owen, this index can be 
computed efficiently by integrating the first derivatives of the reliability function of the system along the main diagonal of 
the unit hypercube. The Samaniego signature of such a system is another important index that can be computed for instance 
by Boland's formula, which requires the knowledge of every value of the associated structure function. We show how the 
signature can be computed more efficiently from the diagonal section of the reliability function via derivatives. We then 
apply our method to the computation of signatures for systems partitioned into disjoint modules with known signatures. 



1. Introduction 

Consider an 71-component system ([n], 0), where [n] = {1, . . . , n} is the set of its components and qb: {0, 1}" -»■ 
{0, 1} is its structure function (which expresses the state of the system in terms of the states of its components). 
We assume that the system is semicoherent, which means that <fi is nondecreasing in each variable and satisfies the 
conditions 0(0, . . . , 0) = and <f>(l, . . . , 1) = 1. We also assume, unless otherwise stated, that the components have 
continuous and i.i.d. lifetimes Ti, . . . , T n . 

Barlow and Proschan |1| introduced in 1975 an index which measures an importance degree for each compo- 
nent. This index is defined by the n-tuple Ibp whose fcth coordinate (fc e [n]) is the probability that the failure of 
component fc causes the system to fail; that is, 

4*> = Pr(T s = T k ) , 

where T$ denotes the system lifetime. For continuous i.i.d. component lifetimes, this index reduces to the Shapley 
value |fT31[T6l , a concept introduced earlier in cooperative game theory. In terms of the values 4>(A) (A £ [n]) of 
the structure function!] the probability 7™ then takes the form 

(1) 4p = E ^tit u{fc})- 0(A)). 

Ac[ n ]\{k} n \\A\) 

The concept of signature, which reveals a strong analogy with that of Barlow-Proschan index above (see Q 
for a recent comparative study), was introduced in 1985 by Samaniego lfT3l[l4l as a useful tool for the analysis of 
theoretical behaviors of systems. The system signature is defined by the n-tuple s whose fcth coordinate Sk is the 
probability that the fcth component failure causes the system to fail. That is, 

s k = Pr(T s = T k:n ) , 

where Tk-. n denotes the fcth smallest lifetime, i.e., the fcth order statistic obtained by rearranging the variables 
Ti, . . . , T n in ascending order of magnitude. 

Boland [3| showed that Sfe can be explicitly written in the form 

(2) -* = E m - t4t e m ■ 

\n-k+l) As[n] \n-k) As[n] 

|-A|=n— fc+1 \A\=n-k 

Thus, just as for the Barlow-Proschan index, the signature does not depend on the distribution of the variables 
Ti, . . . , T n but only on the structure function. 

The computation of lL2 by means of ([T|i may be cumbersome and tedious since it requires the evaluation of 4>(A) 
for every A £ [ti]. To overcome this issue, Owen IfTOlfTTI proposed to compute the right-hand expression in (JTJ only 
from the expression of as a multilinear polynomial function as follows. 
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As a Boolean function, can always be put in the unique multilinear form (i.e., of degree at most one in each 
variable) 

0(x) = £ c(A) Y[xi, 

As[n] isA 

where the link between the coefficients c( A) and the values 4>(A) is given through the conversion formulas (Mobius 
inversion) 

c(A) = £(-l)l A H^(B) and <j>(A) = c(B) . 

B^A B^A 

Owen introduced the multilinear extension of (f> as the multilinear polynomial function 0: [0, 1]™ M. defined by 

0( X ) = y. < a ) n^- 

A£[n] itA 

Example 1. The structure of a system consisting of two components connected in parallel is given by 

4>{X\,X2) = max(si,l2) = X\ II X2 = X\ + X2 ~ X\ X2 , 

where ii is the (associative) coproduct operation defined by xuy = l-(l-x)(l-y). Considering only the multilinear 
expression of <f), we immediately obtain the corresponding multilinear extension <f>(xi , X2 ) = x\ + x-i - x\ x<z. 

In reliability analysis the function <j>, denoted by h, is referred to as the reliability function of the structure </> (see 
12 Chap. 2]; see also lfl2l Section 3.2] for a recent reference). This is due to the fact that, under the i.i.d. assumption, 
we have 



(3) 



F s (t) = ^(t),...,^!)), 



where Fs(t) = Pr(Ts > t) is the reliability of the system and Fk(t) = Pr(Tfc > t) is the reliability of component ft 
at time t. 

We henceforth denote the function cf> by h. Also, for any function / of n variables, we denote its diagonal section 
f(x,...,x) simply by f(x). 

Owen then observed that the fcth coordinate of the Shapley value, and hence the fcth coordinate of the Barlow- 
Proschan index, is also given by 



(4) 



'BP 



[\d k fo{x)dx = f\d k h)( 
Jo Jo 



x) dx . 



That is, Igp is obtained by integrating over [0, 1] the diagonal section of the fcth partial derivative of h. 

Thus, formula (01 provides a simple way to compute I^p from the reliability function h (at least simpler than the 
use of (HJ). 

Example 2. Consider the bridge structure as indicated in Figure Q] The corresponding structure function and its 
reliability function are respectively given by 

<t>(Xl , . . . , X§) = X\ X4 II X2 X§ II X\ X3 £5 II X 2 %3 X4 

and 

h{xi,...,xs) = + + xia:3a;5 + X2a;3a;4 

- X1X2X3X4 - X1X2X3X5 - X1X2X4X5 - X1X3X4X5 - X2X3X4X5 + 2xiX22:3a:4a;5 . 
By using (01 we obtain Ibp = (^)3^)^)M)3^)- Indeed, we have for instance 

4p = f (d 3 h)(x)dx = J X (2x 2 -4x 3 + 2x 4 )dx = ±. 

□ 




Figure 1 . Bridge structure 
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Remark 1. Example[2]illustrates the fact that the reliability function h can be easily obtained from the minimal path 
setfl of the system simply by expanding the coproduct in and simplifying the resulting algebraic expression (using 

Similarly to Owen's method, in this note we provide a way to compute the signature of the system only from the 
reliability function of the structure, thus avoiding Boland's formula (ffl which requires the evaluation of <j)(A) for 
every A £ [n]. 

Specifically, considering the tail signature of the system, that is, the (n + 1) -tuple S = (So, - ■ • , S n ) defined by 
(see ©) 

— 2 1 

(5) S k = £ 8i = — — Z <P( A ) , 

i=k+l \n-k) \A\=n-k 

we prove (see Theorem[5]below) that the coefficient of (x- l) k in the Taylor expansion about x = 1 of the polynomial 

p(x) = x n h(l/x) 

(which is the n-refiected of the univariate polynomial h(x)) is exactly (™) In other terms, we have 

(6) S k = ^-plD k p(l), k = 0,...,n, 
and the signature can be computed by 

(7) s k = S k -i -S k , k = l,...,n. 

Even though such a computation can be easily performed by hand for small n, a computer algebra system can be 
of great assistance for large n. 

Example 3. Consider again the bridge structure as indicated in Figure Q] By identifying the variables xi,...,xs in 
h(xi, . . . , X5), we immediately obtain 

h(x) = 2x 2 + 2x 3 - 5x 4 + 2x 5 , 

from which we can compute 

p(x) = x 5 h(l/x) = 2-5x + 2x 2 + 2x 3 = 1 + 5(x - 1) + 8(x - l) 2 + 2(x - l) 3 , 

or equivalently, 

p(x+l) = l + 5x + 8x 2 + 2a; 3 . 

Using d6) we then easily obtain S = (l, 1, |, |, 0, 0). Indeed, we have for instance (2) S2 = 8 and hence S2 = 4/5. 
Finally, using $7} we obtain s = (0, |, |, |, 0). □ 

This note is organized as follows. In Section 2 we give a proof of our result by first showing a link between the 
reliability function and the tail signature through the so-called Bernstein polynomials. In Section 3 we apply our 
result to the computation of signatures for systems partitioned into disjoint modules with known signatures. 

2. Notation and main results 
Recall that the n + 1 Bernstein polynomials of degree n are defined on the real line by 

&k,»(aO = Qx k (l-x) n - k , k = 0,...,n. 
These polynomials form a basis of the vector space P n of polynomials of degree at most n. 
Proposition 4. We have 

n 

(8) h(x) = Sn-kbk,n(x). 

fc=0 

Thus, the numbers S n -k (fc = 0, . . . , n) are precisely the components of the diagonal section of the reliability function 
h in the basis formed by the Bernstein polynomials of degree n. 



•^Recall that a subset P £ [n] of components is a path set for the function <f> if <j>{P) = 1. A path set P £ [n] is said to be minimal if it does 
not strictly contain another path set. 



Equivalently, (™) Sj. is the coefficient of x k inp(x + 1). 
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Proof. The reliability function can be expressed as 

&(x) = £ 0(A) n«i n (!-«*)■ 

,4c[n] ieA ie[n]\j4 

Its diagonal section is then given by 

= £ 0(A) ^(l-x,)"-^ = t { E 0(A))^(l-x)"- fe 

As[n] fc=0 V |A|=fe / 

and we immediately conclude by (f5]). □ 

By applying the classical transformations between power and Bernstein polynomial forms to Eq. (8), from the 
standard form of h(x), namely h(x) = T,k=o a k % > we immediately obtain 

n-k ( n ~ k \ J \ fc /k\ — 

(9) S t =^«, and Ok = ("JEH) ( l^n-i, fc = 0,...,n 



* (?) 

Remark!. (a) Eqs. ([8]l and © explicitly show that h{x) encodes exactly the signature, no more, no less. 
This means that two ri-component systems having the same h(x) also have the same signature and two 
n-component systems having the same signature also have the same h(x). 

It is also noteworthy that two distinct n-component systems may have the same h(x), and hence the same 
signature. For instance, the 8-component system defined by the structure 

01 (x) = X\ X2 U X2 X3 X4 II X§ Xq X-j Xs 

has the same h(x) as the 8-component system defined by the structure 

02 (x) = xi x 3 u x 2 x± x 5 u x\ x 2 x 6 x 7 x s , 

namely h(x) = x 2 + x 3 - x 6 - x 7 + x s . 

(b) Eq. (O also shows that Sk is the component of h(x) along the basis polynomial b n -k, n - Interestingly, by 
replacing x by 1 - x in (|8), we obtain the following (dual) basis decomposition 

n n 

h(l-x) = Y, S n -kK-k,n(%) = Y S kh,n(x) ■ 
fc=0 fe=0 

(c) Using summation by parts in Eq. dS), we derive the following identity 

n 

fc=l 

where h OSk (x) = SiLn-fc+i ^i,n( x ) is the diagonal section of the reliability function of the (n - k + l)-out- 
of-n system (the structure 0Sfc(x) being the fcth smallest variable £&„). By (0 we see that this identity is 
nothing other than the classical signature-based expression of the system reliability (see, e.g., [14]), that is, 

n 

Pr(T s > t) = Y, s k Pr(T km > t) . 

k=l 

We can now state and prove our main result. Let / be a univariate polynomial of degree m s$ n, 

f(x) = a n x n + ■■■ + ai x + ao . 
The n-reflected polynomial of / is the polynomial f R defined by 

f R (x) = a x n + ai x 71 ' 1 + ■■■ + a n , 

or equivalently, f R (x) = x n f(l/x). 
Theorem 5. We have 

= t(^S k (x~l) k . 

Thus, for every k e {0, . . . , n), the number (?) Sk is precisely the coefficient of (x - l) k of the Taylor expansion 
about x = 1 of the n-reflected diagonal section of the reliability function h. 

Proof. By Proposition!!! we have 

n 

H x ) = Y $ k bn-k,n(x). 
k=0 

The result then follows by reflecting this polynomial. □ 
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From Theorem|5]we immediately derive the following algorithm, which inputs both the number n of components 
and the reliability function h and outputs the signature s of the system. 

Step 1. Express the n-refiected polynomial h R (x) = x n h(l/x) in the basis {(x - l) k : k = 0, . . . , n} 
or, equivalently, the polynomial h R (x + 1) in the basis {x k : k = 0, . . . , n}. That is, 



Cfc x k . 



h R (x) = £c fc (£-l) fc or h R (x + l) = Y, 

k=0 fc=0 

Step 2. Compute the tail signature S : 

Sk = <W(fc) i fc = 0,...,n. 

Step 3. Compute the signature s : 

Sk = Sk-i - S k , k = l,...,n. 

Remark 3. The concept of signature was recently extended to the general non-i.i.d. case (see, e.g., |6l). In fact, 
assuming only that ties have null probability (i.e., Pr(Xi = Tf) =0 for i + j), we can define the probability signature 
of the system as the ?i-tuple p = (pi, . . . ,p n ), where pk = Pr(Ts = Tfc :ji ). This n-tuple may depend on both the 
structure and the distribution of lifetimes. It was proved that in general we have 

n 

(10) Y, Pr(T<j=T i:n ) = £ q(A)<j>(A), 

i=k+l \A\=n-k 

where the function q: 2^ -»■ R, called the relative quality function associated with the system, is defined by q(A) = 
Pr(max ie [„] vA Tj < mm ieA T t ). 

Clearly, the right-hand side of (fTOb coincides with that of © for every semicoherent system when q(A) = 
for every A £ [n] (see [9| for more details). Therefore the algorithm above can be applied to the non-i.i.d. case 
whenever this condition holds, for instance when the lifetimes are exchangeable. 

An n-component semicoherent system is said to be coherent if it has only relevant components, i.e., for every 
k € [n] there exists x e {0, 1}™ such that 0(Ofe,x) + 0(lfe,x), where 0(zfc,x) = tp(x.)\ Xk=z . 

The following proposition gives sufficient conditions on the signature for a semicoherent system to be coherent. 

Proposition 6. Let ([n],</>) be an n-component semicoherent system with continuous i.i.d. component lifetimes. 
Then the following assertions are equivalent. 

(i) The reliability function h is a polynomial of degree n (equivalently, h(x) is a polynomial of degree n). 
(ii) We have 

kodd yft/ t-.— \*/ 



(Hi) We have 



_ , , , , j • 

k odd v a 1 ' k even 

If any of these conditions is satisfied, then the system is coherent. 

Proof. The equivalence (i) o (ii) immediately follows from Theorem[5]and the fact that h(x) is of degree n if and 
only if h R (0) * 0. 

The equivalence (ii) <=>■ (Hi) follows from the straightforward identity 



To see that the system is coherent when condition (i) is satisfied, suppose that component k is irrelevant. Then 
/i(x) = h(lk, x) has less than n variables and therefore cannot be of degree n. □ 

Remark 4. (a) The equivalent conditions in Proposition [6] are not necessary for a semicoherent system to be 
coherent. For instance, the 4-component coherent system defined by the structure 

</>(x) = X\ X2 LI X2 X3 U X3 X4 = X\X2+ X2X3 + X^X^- X1X2X3- X2X3X4 

has a reliability function of degree 3. 
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(b) The 6-component coherent system defined by the structure 

4>i (x) = X\ X2 U X2 %3 X4 II X3 X4 X§ Xq 

has the same h(x) as the 5-component coherent system (or 6-component noncoherent system) defined by 
the structure 

</> 2 (x) = 1113^22:415, 
namely h(x) = x 2 + x 3 - x 5 . We thus retrieve the fact that h(x) does not characterize the system (see 
Remark|2a)) and cannot determine whether or not the system is coherent (see RemarkUa)). 

3. Application: Modular decomposition of system signatures 

We now apply our main result to show that (and how) the signature of a system partitioned into disjoint modules 
can be computed only from the partition structure and the module signatures. 

Suppose that the system is partitioned into r disjoint semicoherent modules (Aj,\j) (j = 1, ■ • ■ ,r), where Aj 
represents the set of the components in module j and \y {0, l} Aj -»■ {0, 1} is the corresponding structure function. 
Let rij denote the number of components in Aj (hence Yij=i n j = n ) ar, d l et Sj = (Sj.o, . . . , Sj. nj ) denote the tail 
signature of module j. 

If V> : {0, l} r -*■ {0, 1} is the structure function of the partition of the system into modules, the modular decompo- 
sition of the structure <f> of the system expresses through the composition 

<Kx) = ^(xi(x Al ),...,x,(x A 0), 

where x Aj = {xi)^A (see 12 Chap. 1]). Since the modules are disjoint, this composition extends to the reliability 
functions h$, h^, and h Xj of the structures 0, ijj, and Xj> respectively; that is, 

(11) Mx) = K{Ki(* Al )>•••, h Xr (* Ar ))- 

Indeed, the right-hand side of ( fTTT i contains no powers and hence is a multilinear polynomial. 

According to Theorem|5] the tail signature of the system can be computed directly from the function 

(12) h$(x) = x n h^(x-^h^(x),...,x- n ^hl(x)), 
where h?: is the n, . -reflected of the diagonal section of h v , that is, 

(B) h*{x) = x^' l ;)s J . l ,(.r \f. 

Interestingly, Eqs. dT2b and (TT~3b show that (and how) the signature of the system can be computed only from the 
structure t/j and the signature of every module. Thus, the complete knowledge of the structures \i, ■ ■ ■ , Xr is n °t 
needed in the computation of the signature of the whole system. 

Example 7. Consider a 7-component system consisting of two serially connected modules (hence ip(zi, z%) = z\Z-x) 
with signatures Si = (~, |, 0J and s 2 = (0, |, |, 0), respectively. By ( TT3l we have 

(x) = 2x-l and /j* (x) = 2x 2 - 1 . 

By d 1 2b we then obtain 

hf(x) = x 7 (x- 3 (2x- l)x- i (2x 2 - 1)) = l-2x-2x 2 + 4x 3 , 

from which we derive the system signature s = , 0, 0, 0). □ 

As an immediate consequence of our analysis we retrieve the fact (already observed in (8]; see |@][5) for earlier 
references) that the signature always decomposes through modular partitions We state this property as follows. 

Theorem 8. The signature of a system partitioned into disjoint modules does not change when one modifies the 
modules without changing their signatures. 

A recurrent system is a system partitioned into identical modules. Thus, for any recurrent system we have ri\ = 
•■■ = n r = n/r and \i = = Xr = X- Eq- (TT~2b then reduces to 

(14) h$(x) = x n h^(x- n l r h*{x)) = h*(x) r h*(x n ' r h^(x)- 1 ) . 

Thus, to compute the tail signature S of the whole system the knowledge of the structures i/j and \ can t> e simply 
replaced by the knowledge of their corresponding tail signatures S^, and S x , respectively. 

**This feature reveals an analogy for instance with the barycentric property of mean values: The arithmetic (geometric, harmonic, etc.) mean 
of n real numbers does not change when one modifies some of the numbers without changing their arithmetic (geometric, harmonic, etc.) mean. 
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Example 9. Consider a system partitioned into r modules, each of whose consists of two components connected in 
parallel (system with componentwise redundancy). In this case we have h^(x) = 2x - 1 and 

By (O it follows that 
from which we derive 

[t/2\ (r\(r-k\ 

"c V \k \t-2k n £-2k ~c n s s n 

be = l_, — ?27\ — 2 b -4>,k, 0^£^2r. 

fc=max(<?-r,0) \ £ ) 
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